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I .  INTRODUCTION 

A  recent  test  involving  an  air  defense  missile  system  required  the  con¬ 
sideration  of  low  altitude  target  detectability  and  clutter  and  multipath 
effects  for  both  the  ground-based  fire  control  radar  and  the  missile  RF  seeker. 
A  visit  to  the  test  site  at  White  Sands  Missile  Range  (WSMR),  New  Mexico,  and 
discussions  with  personnel  familiar  with  the  area  suggested  that  the  targets 
might  indeed  be  flown  below  the  radar  horizon  in  some  regions  of  interest,  and 
that  portions  of  the  terrain,  masked  from  the  radar  or  seeker,  might  change 
clutter  and  multipath  predictions. 

An  initial  analysis  was  done  using  contour  maps  of  the  region,  followed 
by  a  visit  to  the  site.  These  two  approaches  did  not  exactly  agree,  since  the 
on-site  examination  indicated  a  more  severe  problem  than  predicted  from  the 
contour  maps.  This  may  have  been  due  to  difficulty  in  comparing  the  contour 
data  representing  fairly  large  intervals  to  the  visual  impact  of  the  site. 

Some  difficulty  was  noted,  however,  in  comparing  the  map  features  with  the 
actual  terrain.  In  any  event,  a  more  exact  method  was  required  and  the  avail¬ 
ability  of  digitized  terrain  elevation  data  suggested  a  computer  masking 
analysis. 

A  set  of  computer  programs  was  implemented  to  reformat  the  data  and  per¬ 
form  a  line-of -sight  masking  analysis.  This  document  presents  the  results  cf 
a  radar  masking  investigation  of  two  separate  locations  at  WSMR.  Described  in 
the  following  sections  are  the  terrain  type  and  location,  and  data  analysis 
and  manipulation.  Results  are  presented  in  the  form  of  ground  masking  plots. 

A  sample  computer  program  is  included  in  the  Appendix. 


II.  DESCRIPTION  OF  TERRAIN 

The  testing  area  is  located  within  the  White  Sands  Missile  Range,  New 
Mexico.  The  two  radar  locations  considered  are  marked  EAST  CENTER  50  (EC  50) 
and  WEST  CENTER  50  (WC  50)  on  the  map  of  Figure  1 .  The  sector  of  interest  was 
roughly  rectangular  in  shape  and  bounded  by  106  16'  and  106  26'  longitude,  and 
33°8'  and  33°25'  latitude.  These  limits  are  indicated  in  Figure  1.  Target 
flights  were  to  be,  in  general,  along  north-south  paths  within  this  region. 

The  locations  of  the  two  sites,  as  used  in  the  analysis,  are  given  in  Table  1. 


Table  1.  Radar  Locations 


Location 

Longitude 

Latitude 

EC  50 

WC  50 

106o17'3r' 

1 06°26 ' 1 1 . 5" 

33°8  *21" 

33°8'43" 

The  surface  characteristics  of  the  test  site  region  were  described  in 
three  forms: 

•  Visual  from  a  visit  to  the  site 
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•  Geological  survey  topographic  maps 

•  Defense  Mapping  Agency  (DMA)  terrain  altitudes  on  magnetic  tape. 

Each  data  source  aided  in  estimating  the  radar  coverage. 

Gross  terrain  features  are  indicated  in  Figure  1.  The  terrain  for  both 
sites  is  gently  rolling  to  flat  within  most  of  the  test  region.  A  Malpais 
region  exists  to  the  north,  with  the  San  Andres  Mountains  to  the  northwest. 

A  dry  stream  bed  runs  north  to  south  over  part  of  the  western  edge  of  the  area. 

A  gentle  rise  was  visually  noted  at  5  or  6  miles  from  the  radar  site  at  EC  50 
and  appeared  to  mask  much  of  the  terrain  in  the  test  area.  This  was  not  evi¬ 
dent  on  the  contour  map  and  led  to  concern  over  the  extent  of  masking  and  other 
possible  "hidden"  regions.  These  questions  led  to  the  analysis  of  DMA  digitized 
terrain  data  on  a  CYBER  74  digital  computer. 


III.  MASKING  ANALYSIS 

The  DMA  terrain  data  consisted  of  terrain  elevation  above  mean  sea  level 
(MSL)  truncated  to  the  nearest  meter.  The  data  files  are  arranged  into  I- 
degree  by  1 -degree  geographic  areas,  with  each  file  containing  samples  falling 
into  a  single  1-degree  square.  The  reference  for  each  file  is  the  southwest 
comer  of  the  degree  square.  Each  file  contains  1200  data  records,  each  having 
a  constant  longitude  value.  The  first  value  in  each  record  is  the  southern¬ 
most  elevation  within  the  square;  the  last  value  is  the  northernmost  elevation. 
Each  such  record  contains  1200  values  so  that  the  terrain  elevation  is  sampled 
on  3-second  intervals  in  longitude  and  latitude.  On  an  earth  model  with  a 
4,000-foot  MSL  elevation,  this  sample  spacing  is  approximately  250  feet  in 
longitude  and  300  feet  in  latitude  for  the  region  of  interest.  These  data 
require  several  transformations  to  be  useful  for  a  radar  terrain  analysis. 

Consider  the  radar  radial  elevation  cut  of  Figure  2.  Transformation  of 
the  terrain  elevation  to  radial  elevation  cuts  allows  a  rather  simple  algorithm 
to  search  for  the  masked  regions,  shown  as  shaded  areas,  and  to  store  the  start 
and  stop  ranges  of  masking.  Repeated  radial  profiles  then  produce  a  map  simi¬ 
lar  to  Figure  3  where  the  shaded  regions  now  represent  the  surface  area  illum¬ 
inated  by  the  radar.  The  negative  of  this  presentation  could  also  be  used  to 
represent  the  masked  area. 

The  data  transformation  and  manjnulations  required  are  listed  below: 

•  Transformation  from  geodetic  coordinates  (latitude,  longitude,  eleva¬ 
tion)  to  geocentric  coordinates  (Xc,  Yc,  Zc) . 

•  Transformation  from  geocentric  coordinates  (Xc,  Yc,  Zc)  to  tangent 
plane  coordinates  through  a  point  on  or  above  the  surface. 

•  Reverse  order  the  data  left  to  right  (EC  50  only). 

•  Interpolate  in  three  dimensions  to  produce  radial  elevation  profiles. 

•  Masking  algorithm  along  radial  elevation  profiles  to  lv cate  masked 
regions . 


Figure  2.  Radar  masking  definitions. 


«  Graphic  presentation  of  results. 

The  two  transformation  algorithms  were  taken  from  Reference  ClD  and 
applied  directly  to  the  data.  No  difficulties  were  encountered;  however,  it 
should  be  noted  that  the  Cartesian  systems  are  "left-handed."  For  example,  in 
the  tangent  plane  system,  the  X-axis  lies  along  a  north-south  line  through  the 
tangent  point  and  is  positive  north  of  the  point;  Z  is  altitude  normal  to  the 
tangent  plane  and  positive  up;  and  Y  is  positive  to  the  east.  The  ellipsoidal 
earth  model  used  was  the  Clarke  Spheroid  of  1866  with  semi-major  axis 
20,925,832,136’  and  eccentricity  (e)  0.0822718536. 

Item  3,  reversal  of  the  data,  is  required  for  EC  50  due  to  the  order  of 
the  data.  The  DMA  data  start  in  the  southwest  corner  of  the  region,  hence  the 
regions  at  longest  range  from  the  radar  appear  first  on  the  tape.  The  data 
were  reordered  such  that  near  ranges  were  physically  located  first  on  the  data 
file.  The  interpolation  routine  could  then  read  a  pair  of  records,  perform 
the  required  interpolation  to  convert  X,  Y,  Z  to  range,  azimuth,  elevation, 
store  these  values,  and  read  in  the  next  set  of  elevation  data.  The  mask 
algorithm  was  then  used  to  search  along  constant  azimuth  profiles  for  hidden 
elevation  regions  and  the  results  plotted  with  routines  developed  locally  and 
reported  in  Reference  [2]. 

The  WC  50  data  did  not  require  reversal,  so  Item  3  was  not  used.  The 
remainder  of  the  analysis  was  identical. 


IV.  RESULTS 

The  graphic  presentations  of  radar  masking  resulting  from  this  analysis 
are  given  in  two  forms:  the  range-azimuth  "map"  described  earlier,  plus  a 
series  of  elevation  profiles  to  aid  in  interpretation  of  the  map  and  estimation 
of  masked  flight  altitudes  for  the  targets.  Three  radar  altitudes  of  interest 
are  also  presented;  each  represents  the  height  of  the  radar  antenna  center 
above  the  local  terrain.  Figures  4,  5,  and  6  are  maps  for  EC  50  of  illuminated 
regions  for  the  three  heights  for  the  azimuth  coverage  of  -5°  to  -70°.  Range 


8 


CROSS  RANGE 


Figure  3-  Masking  Diagram 


coverage  was  0  to  100,000  feet,  with  the  cross  range  dimension  limited  to 
70,000  feet  for  this  presentation.  Azimuth  increments  were  1  degree;  range 
increments  50  feet. 

As  expected,  reduction  of  antenna  height  significantly  reduces  the  illum¬ 
inated  terrain  area.  Clutter  would  be  reduced  somewhat  for  medium  ranges,  but 
the  mountain  ranges  at  long  range  might  still  produce  significant  clutter. 

Radial  elevation  profiles  for  EC  50  are  presented  in  Figures  7  through  19 
for  the  40-foot  radar  altitude  only.  Other  altitudes  may  be  investigated  by 


■L1  ,.l".lUI 


shifting  the  left  hand  altitude  scale.  The  mountains  appear  quite  abrupt  due 
to  the  difference  in  scale  factors  on  the  two  axes.  Also,  note  the  "tilt"  of 

O  o 

flat  areas  due  to  earth  curvature.  In  the  region  -20  to  -40  a  target  would 
have  to  be  below  75  feet  to  100  feet  altitude  for  a  40-foot  radar  altitude  in 
order  to  be  masked,  and  then  only  at  ranges  of  55  to  65  kilofeet;  however,  for 
a  9-foot  radar  altitude,  the  same  target  would  be  masked  at  about  125  feet  alti¬ 
tude.  These  boundary  altitudes  between  masked  and  unmasked  could  be  signifi¬ 
cantly  higher  when  vegetation  is  added  to  the  terrain  and  other  propagation 
problems  are  considered. 

A  similar  set  of  data  for  WC  50  is  given  in  Figures  20  through  36.  For 
this  location,  only  the  9-foot  antenna  height  was  used,  and  an  investigation  of 
clutter  fences  was  done  for  this  site.  Figures  20  and  21  are  illuminated 
regions  in  two  scale  factors  without  the  fence,  while  Figures  22  and  23  show 
the  effects  of  a  circular  section  clutter  fence  whose  location  was  100  feet 
from  the  radar  and  whose  top  was  the  same  as  the  antenna  feed,  9  feet.  The 
reduction  of  close-in  illuminated  area  is  significant,  without  noticed  increase 
in  long-range  masking.  Figures  24  through  36  are  radial  altitude  profile  cuts 
similar  to  those  for  EC  50. 

Note  that  in  Figures  20  through  23,  the  algebraic  sign  of  the  cross  range 
axis  has  been  changed  to  facilitate  plotting.  The  sign  conventions  of  section  3 
were  used  in  developing  the  figure. 

The  possibility  of  terrain  masking  of  low  altitude  targets  is  lower  for 
this  site.  For  most  of  the  sector,  the  long  range  clutter  problem  has  been 
significantly  reduced  since  the  long  range  mountain  slopes  only  appear  almost 
due  north  from  the  site.  In  this  sense,  WC  50  presents  a  more  typical  clutter 
environment  than  EC  50. 

Some  caution  should  be  exercised  in  using  the  results,  in  particular  the 
clutter  fence  data.  The  computer  algorithm  claculated  line-of -sight  masking  in 
which  the  edge  diffraction  effects  were  ignored.  Therefore,  to  a  certain  extent, 
many  of  the  shadow  regions  will  be  partially  illuminated.  This  is  believed  to 
be  an  insignificant  effect  for  the  terrain  masking,  since  many  other  factors, 
such  as  surface  roughness  and  vegetation,  will  reduce  the  abruptness  of  the  edge. 
This  is  not  the  case  for  the  conducting  clutter  fence.  A  recent  analysis  of 
this  problem  [3]  indicates  that  no  more  than  a  few  dB's  reduction  of  the  clutter 
energy  should  be  expected  for  this  geometry.  Therefore,  the  clutter  fence  was 
not  constructed  and  no  evaluation  is  possible. 


V.  CONCLUSIONS 

The  slight  rise  evident  from  visual  inspection  was  quite  obvious  on  the 
DMA  data  and  did  mask  much  of  the  region  from  20  kilofeet  range  out  to  the 
mountains  at  60  to  80  kilofeet  for  EC  50.  This  feature  is  not  that  evident  on 
the  contour  maps,  so  that  the  computer  analysis  did  predict  some  of  the  masking 
evident  in  actual  flights.  Target  altitudes  for  masking  did  not  agree  as  well 
as  expected,  but  vegetation  effects  were  not  included  in  this  analysis  and  the 
precision  of  reported  target  altitude  is  unknown. 
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Figure  4.  Illuminated  region  from  EC  50,  40-foot  radar  height. 
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Figure  9.  EC  50  elevation  prof 
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Figure  12.  EC  50  elevation  profile,  -35  azimuth,  40-foot  radar  height. 


Figure  14.  EC  50  elevation  profile,  -45  azimuth,  40-foot  radar  height 


Figure  15.  EC  50  elevation  profile,  -50  az 
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Figure  16.  EC  50  elevation  profile,  -55°  azimuth,  40-foot  radar  height 
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Figure  17.  EC  50  elevation  profile,  -60  azimuth,  40-foot  radar  height 
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Illuminated  region  from  WC  50,  expanded  scale,  9-foot  radar  height. 
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Figure  24.  WC  50  elevation  profile,  5  azimuth,  9-foot  radar  height. 


Figure  25.  WC  50  elevation  profile,  10  az 


Figure  26.  WC  50  elevation  profile,  15  azimuth,  9-foot  radar  height. 
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Figure  27.  WC  50  elevation  profile,  20  azimuth,  9-foot  radar  height 
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Figure  29.  WC  50  elevation  profile,  30°  azimuth,  9-foot  radar  height 


fC  50  elevation  profile,  35  azimuth,  9-foot  radar  height 


Figure  31.  WC  30  elevation  profile,  AO  azimuth,  9-foot  radar  height. 
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APPENDIX 

COMPUTER  PROGRAMS 


47 


Several  computer  programs  were  developed  in  the  process  of  evaluating 
approaches  and  performing  the  analysis.  In  general,  the  tasks  of  transfor¬ 
mation,  interpolation,  and  masking  are  the  significant  steps.  Since  the 
transformation  programs  are  documented  in  the  reference  they  will  not  be 
repeated  here.  The  only  note  of  caution  to  the  prospective  user  is  that  the 
sense  of  the  various  coordinate  systems  must  be  clearly  understood. 

The  interpolation  and  masking  steps  will  be  represented  by  an  example 
program.  This  routine  used  as  input  on  Tape  10  data  already  transformed  to 
radar  tangent  plane.  The  program  selects  a  radial  for  analysis  and  inter¬ 
polates  data  points  in  50-foot  increments  via  subroutine  TINTER.  After 
creation  of  the  radial  elevation  profile,  subroutine  MASKK  is  called  to  do 
the  masking  calculations.  Start  and  end  points  of  illumination  are  con¬ 
nected  by  a  solid  line  and  plotted  using  subroutine  PLTTEK.  This  utility 
to  produce  CRT  displays  of  the  data  is  also  documented  in  the  references 
and  not  included  here. 

Outputs  from  this  program  include  the  CRT  display  of  illuminated 
regions  with  the  same  data  stored  on  magnetic  tape  (TAPE  2)  and  signifi¬ 
cant  comments  concerning  the  results  on  OUTPUT. 

This  description  is  intended  to  give  the  flavor  of  the  computations 
rather  than  an  exhaustive  description.  No  extensive  attempt  was  made  to 
maximize  the  computational  speed  or  efficiency  of  these  programs.  If  their 
use  is  contemplated  in  a  production  environment,  such  improvements  should 
definitely  be  considered. 


PPOGRAM  MASK 


If 


IS 


20 


2S 


30 


35 


40 


45 


SO 
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PROGRAM  MASK  <  INPUT  .01 ITPI  IT  .  T  APFS=  INPUT  .  T APE4=0UTPUT  . 

.  TAPF10.TAPE1.TAPE2) 

f 

r 

DIMENSION  IHEADUO)  .XA<340)  .Y  A  <  340 )  .7M  340)  .X8(340)  .YB1340)  . 
.  ZB(340) ,H(2000) .o(?000) 

DIMENSION  LZ<2> 

DATA  PI/1.14159/ 

r 

r 

WRITE  11 

11  FORMAT ( 1H1 ) 

READ ( 10) IMF AD 
WRITE (?)IHEAD 
WRITE  l.IMEAO 

1  FORMAT (1X.10A10) 

r 

c  angle  in  degrees 

c 

READ  2. THRA Y .HO .CUTS 

2  FORMAT (3F10.2) 

WRITF  2.THRAY «H0 .CUTS 
ICUTS=CUTS*.5 
WRITE (2) THRA Y. HO. CUTS 
C 
c 

DO  1  OR  1  =  1.2000 
1  OR  R(I)  =  (I-1)*50. 

L  Z (2)  =  10H* 

THRAYOrTHRAY 
N340=340 
N2000=2000 
lx=iohrangf  ft* 

LY=10HALT  TP* 

c 

r 

XMN=0. 

YMN=YMX 
XMX=  100000. 

DO  500  IRAY=1  .  I  Cl  IT  S 
00  12  1=2.2000 

12  H(I)=0. 

THRAY=THRAVO* ( IRAY-1 )*1 . 

REWIND  10 

READ!  10) IHEAD 

THRAYD=THRAY 

THRAY*THRAY*PI/150. 

TANTH*TANITHRAY) 

COSTHRaCOSITHRAY) 

SINThR=SIN(THRAY) 

H  ( 1  )  *0  . 

P(1)*0. 

WRITE  10.IRAY.THRAYD 

10  FORMAT!/, IX. *1RAY  =  *.I4,*  THRAYD  a  *.F10.2) 

DO  3  1=1.200 
I  Y=  I 

READ ( 1 0) XA.YA.7A 


000410 

000420 

000430 

000440 

000450 

000440 

000470 

000450 

0004QO 

000500 

000510 

000520 

000530 

000540 

000550 

000540 

000570 

OOOSRO 

O005R0 

000600 

000610 

000620 

000630 

000640 

000650 

000640 

000670 

000650 

000690 

000700 

000710 

0007)5 

000716 

000717 

000720 

000730 

000740 

0007S0 

000740 

000750 

000755 

000756 

000790 

000500 

OOOSIO 

000820 

000830 

000840 

000850 

000840 

000870 

000850 

000855 

000856 

000890 

000900 

000910 
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program  mask 


IFtYA(l) .LT.O.)GO  TO  4 
1  CONTINUE 

60  4  CONTINUE 

BACKSPACE  10 
IY=IY-1 
IX  =  I 

YAM)  =  Y8  ( 1 ) 

6S  WRITE  60«IY 

60  FOPMAT!*  7FRO  COLUMN  r  *.IS) 

DO  110  1=2.2000 

px=<i-n*so. 

XD=RX*COSTHR 
70  YD=PX*SINThR 

IF  (Y0.GT.Y8!  IX)  .  AND.  Y!).(_T.  VA  !  I  X  >  )  C>0  TO  10? 

DO  S  J=1 .200 
JY=J 

REAOM0IXA.YA.7A 

75  IF(EOFMO)  11000.100 

100  CONTINUE 

1F(YAMX).LT.Y0)C,0  TO  6 

5  CONTINUE 

6  CONTINUE 

flo  IY=IY»JY 

BACKSPACE  10 
BACKSPACE  10 
PEAOMO)  XA.YA.ZA 
IF (EOF  <10)1 1000. 101 
AS  101  CONTINUE 

PEAOMO)  XB.Y8.7B 
IF (EOF (10) 11000.102 
102  CONTINUE 

IF  ( I  «LT  .2)  WPITF  M.IY.I 

90  6j  FORMAT!*  COLUMN  *.IS.*  ChOSFN  FOP  PASS  *.IS) 

00  7  J=1,1S0 
JX  =  J 

IFIXAU)  .GT.XOIGO  TO  8 

7  CONTINUE 

PS  A  CONTINUE 

IX=JX-1 

IYY=IY-1 

CALL  T  INTER  (XA. YA. 7A. XB. YB, ZB. IX. IYY.XD.YO.HM)  .N140.I) 
110  CONTINUE 

100  GO  TO  1002 

1000  CONTINUE 
WRITF  1001,1 

1001  FORMAT!*  FOF  ENCOUNTERED  TN  PASS  *.IS) 

1002  CONTINUE 

10S  WRITE(2)TMRAYD.  N2000.R.H 

CALL  MASKK (P.H.2000.THPAY, IPAY.XMN.XMX.YMN.YMX) 

ENCODE M0.90.LZM)  1THRAYD 
90  FORMATIF10.2) 

CALL  PLTTEK1R.M.2000.0.LX.LY.L7.XMN.XMX.YMN.YMX) 

110  SOO  CONTINUE 

65  FORMAT !  RE  15.7 1 

ISO  FORMAT IFIO. 2) 

CALL  FINITT 
STOP 

IIS  END 


000920 

000910 

000940 

0009S0 

0009*0 

000970 

0009»0 

000990 

001000 

001010 

001020 

001010 

001040 

ooioso 
001060 
001070 
0010A0 
001090 
001100 
001110 
001120 
001110 
001140 
0011S0 
0011*0 
001170 
0011A0 
001190 
001200 
001210 
001220 
001210 
001240 
001290 
001260 
001270 
001 2A0 
001290 
001300 
001310 
001320 
001310 
001340 
0013SO 
001360 
001370 
O013A0 
001390 
001400 
001401 
001402 
001403 
001410 
001440 
001490 
001520 
001570 
001540 
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SURPDUTINE  MASKK 


l 


5 


in 


15 


20 


25 


30 


1* 


40 


SUBROUTINE  MASKK (R.H.NN. TMruT . TP . XMN. XMx . Y“N. YMX 1 

0032)0 

DIMENSION  R(NN) ,H(NN).X(2> .Y(2) 

003220 

NNS=0 

00323" 

cccc 

IE(IP.GT.l>GO  TO  BO 

003240 

GO  TO  BO 

003241 

TMXX=AflS<YMXI 

0032*0 

X<1)=XMN 

003240 

Y(1)=YMN 

003270 

LX=10HO«N  PNG* 

0032«0 

LY=10HCROSS  PNG* 

0032B0 

LZMONMASK* 

003300 

CALL  PLTTEK (X.Y, 1 .NNS.LX.Lv.L7.XMN.XMX, YMN. Y»XX) 

0033)0 

<30 

CONTINUE 

003320 

c 

003330 

R1=RU) 

003340 

I  SET =0 

003350 

TANTHX=-100. 

003340 

00  100  1*2, NN 

003370 

TANTM=MU)/R(I) 

0033B0 

IFITANTH.LT. TANTHX1GO  TO  I'd 

0033B0 

rrasET.Eo.i»Pi*pf  n 

003400 

TANTHX=TANTH 

003410 

ISFT=0 

003420 

GO  TO  100 

003430 

101 

CONTINUE 

003440 

IE ( I SET.EO.O ) PM1=P  < I ) 

0034*0 

IEdSET.EO.llGO  TO  100 

003440 

ISET=1 

503470 

R2=R(I-I) 

0034BO 

X(1 )=P1*C0S(THCUT) 

0034RO 

Y(1)=R1»ABS(SIN<TnCUT) ) 

003500 

X  (2)  =R2*C0S  (THCllTl 

003510 

Y(2)=R2*SIN(THCUT) 

003520 

Y (?) =ABS ( Y 12) ) 

003530 

IE(R1.F0.R?)NRITF  40 ,R1 .R2. IP. I . ThCUT 

003540 

cccc 

CALL  PLTTCK(X.y.~S.mS> 

003550 

100 

CONTINUE 

003540 

40 

FORMAT (*  ERROR  IN  M»SK  *.?E1 ?.?.?|«,Ei o.?) 

003570 

RETURN 

003500 

END 

003SB0 
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SURPOUTINE  TINTER 


1 


s 


IP 


is 


?<> 


2S 


10 


SUBROUTINE  TINTERIXA.YA.ZA.XR.YR.ZR.IX. ir.XD.YO.77.NN. t ) 
DIMENSION  XA (NN) .YA (NN>  «7A INN) .XR !NN) .YR  <NN) .ZR INN) 
0ELX=X4<IX.1>-XA(IX| 

0ELY1=YR(TX)-YA(IX) 

DEL  Y2=YRdX.ll-YA < IX*1 > 

DYl=Yn-YA(IX) 

DY?sYO-YA(IX*l) 

AA1=!ZR(IX)-ZA(IX)  )«nYt/0£l  Y|.7AdX> 

AA2=(Z5<IX.1>-ZA<IX*1)  )*0Y2/DFLY2.ZA<IX.l) 

OX=XO-XA(IX) 

ZZ=(AA?-AA1)»0X/0ELX*AA1 

IE (Z7.LT. -1000. I  WRITE  ?, IX. IY.XA (IX) .XA  » I X* 1 ) .XR(TX) .XR ( I  X ♦ 1 ) . 
.  YA(IX).YAdX*l>.YSdX).YR(!X.l). 

.  ZA(1X>  .ZAdx*l)  .ZR(IX)  .ZRdXd)  . 

.  OELX.OELYl.OELY?.OY| ,ny?.AAI .AA2.DX.Z? 

?  F0RMATdX.2IS.12E5. I ./ . I x .oE)2.4 ) 

ir(xn.LT.xAdx)  .or.xo.ot.xa  (ix»i )  > write  io.ix.iy. 

.  AA1 .AA2.XA (IX) .XA (IX.l) .XD.I 

10  FORMAT (*  FIRST  ERROR  *.  2IS.SE1S.7.*  PASS  *.I5> 

IF  (XO.LT.XB  (IX)  ,OR.XO.ST.XA(IX*1 )  )  WRITE  U.IX.IY. 

.  AA1.AA2.XR (I») ,XR( I  X  *  1 ) ,XD,I 

11  FORMAT (•  SECOND  ERROR  *.2  IS. SE IS. 7.*  PASS  *.I5> 

IF(YO.OT.YAdX)  .OR.YO.LT.YR(IX)  )  WRITE  12.IX.IY.AA)  .AA?. 

.  YAdX)  .YR<IX)  .YO.I 

12  FORMAT!*  TH I RO  ERROR  *.2 1 S ,SF IS . 7 , •  PASS  *.IS) 

IF  (YO.fiT.YAdX.il  .OR.YO.LT.YR(IX*I  )  )  WRITE  13.IX.IY.AA]  .A  A?  . 

.  YA(IX.  1)  .YRdX.l)  .YO.I 

13  FORMAT!*  FORTH  ERROR  *.2TS,SF1S.7.*  PASS  *.151 

RETURN 

END 


001SS0 
001550 
001570 
0015R0 
0015RP 
001600 
091610 
001620 
001630 
001640 
001650 
001651 
0016S? 
001653 
001654 
001655 
001640 
001670 
0016R0 
0016R0 
001700 
001710 
001720 
001730 
001740 
001750 
001740 
001770 
001 7»0 
0017R0 
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